Atomic current in optical lattices: Esaki-Tsu equation revisited 
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The paper discusses the master equation approach to derivation of the Esaki-Tsu equation for 
drift current. It is shown that the relaxation term in the master equation can be identified by 
measuring the velocity distribution of the carriers. We also show that the standard form of the 
relaxation term, used earlier to derive Esaki-Tsu equation, predicts unphysical velocity distribution 
and suggest a more elaborated relaxation term, which is argued to correctly capture the effect of 
bosonic bath in experiments on atomic current in optical lattices. 



1. Recently much attention has been payed to the 
transport of cold atoms in optical lattices [J H, S 0, H, 
S S S| ■ This is not only due to the fact that this system 
mimics electrons in a crystal lattice, thus allowing stud- 
ies a variety of transport phenomena of the solid state 
physics. More importantly, the system of cold atoms in 
optical lattice offers a possibility for direct measuring of 
the quantities, which were earlier measured only indi- 
rectly, and opens an access to the parameter regimes, 
which were inaccessible with solid crystals. This ad- 
dresses a number of questions, which may require a re- 
vision of the results obtained earlier in the field of the 
condensed matter physics. In this brief report we discuss 
the Esaki-Tsu equation for the drift current, suggested by 
Esaki and Tsu in 1970 with respect to the biased semi- 
conductor super-lattices 0] . Recent experimental 1] and 
theoretical 0, Q studies of the atomic dynamics in tilted 
optical lattices indicate that this equation also holds, at 
least qualitatively, for the cold atoms. The aim of the 
present work is to explore the 'degree of uncertainty' in 
the Esaki-Tsu equation which, as it will be shown below, 
comes from different models for interactions between the 
system (the carriers, in what follow) and its environment 
(the bath). We shall also show that in the case of cold 
atoms one can get information about the actual interac- 
tions by measuring the velocity distribution of the carri- 
ers instead of measuring the net current. 

2. First we recall the reader the quantum-mechanical 
derivation of the Esaki-Tsu equation The starting 
point is the master or Liouville equation on the reduced 
density matrix of the carriers. 
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where H is the carriers single-particle Hamiltonian in the 
tight-binding approximation. 
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and C{p) the relaxation term, which takes into account 
the effect of a bath. (To simplify equations we set the 
lattice period and Planck's constant to unity.) It is fur- 
ther assumed that that for F = the bath brings system 



into the thermal state po and that this process is charac- 
terized by a single relaxation constant 7, i.e.. 



^(p) = -lip - Po) , Po ^ exp(-/3i?o) 
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Clearly the thermal density matrix po is diagonal in the 
quasi-momentum basis |fc) ~ L^^/^ exp(i27rfc^/L)|Z). 

On the other hand, the Hamiltonian H is diago- 
nal in the basis of the Wannier-Stark states |to) = 
J2i Ji-miJ I F)\l) (here Jn{z) are Bessel functions of the 
first kind), with eigenvalues forming the Wannier-Stark 
ladder. Using the latter basis the stationary solution 
of the master equation ([1]) with the relaxation term ([3|) 
reads [ill. 
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where I„(z) are the second kind Bessel functions. Fi- 
nally, substituting Q into expression for the mean car- 
rier velocity, 

v = Tt[vp], v = lj2i\l + l){l\-h.c.) , 
I 

one recovers the Esaki-Tsu equation with the tempera- 
ture dependence given in the prefactor f{(3), 
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At this point we note that the result ([5]) heavily relies 
on a particular choice of the relaxation term ([3]), which 
by no means can be considered as well justified. For this 
reason we consider a more general relaxation term, 



(6) 



which is explicitly written here in the Lindblad form [12| . 
In Eq. ([6]) the operators the quasimomentum 

representation are given by the matrices 
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and the coefRcients W{s, q) have sense of the transition 
rates between different quasimomentum states in the ab- 
sence of static forcing. The only fundamental restriction 
on the coefficients W{s, q) is 



Wiq,s)=eP^^^-^'^Wis,q) , Eg 
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(8) 

which insures relaxation of an arbitrary initial state to 
the thermal equilibrium. In all other aspects the entries 
W{s, q) are model specific. In the rest of the paper we 
analyze two of these models, which are especially useful 
for understanding physics behind Eq. ([5]). 
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3. It is convenient to rewrite the master equation ([T]) 
with the relaxation term ^ as an equation on matrix 
elements of the density matrix in the quasimomentum 
representation, 
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First we consider the case where the transition rates are 
given by W{s,q) = dg.o- Then, assuming for a moment 
F = 0, the balance equation (fTO| simplifies to 
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Physically this corresponds to zero temperature of the 
bath and selection rules, where any quasimomentum 
state relax directly into the ground state. We would 
like to stress that the considered choice of the coefficient 
VF(s, q) is the only case where the relaxation term ^ 
coincides with ([3]). For all other choices a reduction of 
(|6l) to the simple form ([3]) requires some approximations 
or is impossible in principle. 

Let now F ^ Q. Assuming the limit of an infinite lat- 
tice, where pk.k{i) — p{k,t) is a continuous function of 
the quasimomentum, the coherent evolution term modi- 
fies Eq. (HI]) as 

(12) 

Since J p{k, t)dk ~ 1 due to the normalization, the sta- 
tionary solution of (fT2|) is given by 
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with discontinuous jump at fc = 0. (Note that p{k,t) is a 
periodic function of the quasimomentum.) Substituting 



FIG. 1: Numerical analysis of the model No.l. The mean 
velocity as the function of time (solid line in the upper panel) 
and stationary distribution over the quasimomentum (asterisk 
in the lower panel), compared with the analytical result (dash- 
dotted line). Parameters are F = 0.2, 7 = 0.04, ksT = 0, and 
L = 64 (periodic boundary conditions). For the sake of com- 
parison dashed lines show the mean velocity and stationary 
distribution at infinite temperature. 



(10) this stationary solution into expression for the mean ve- 
locity, V = J^j. sin(27rfc/L)pfe^fe, we recover Eq. ([5]) for 
the net current. 

As an illustration to the above analysis and to test the 
code Fig. [1] and Fig. [5] show numerical solution of the 
master equation for VF(s, g) — Sq^s- In our numerical ap- 
proach we solve the equation in the Wannier basis for 
a finite lattice with periodic boundary conditions. Note 
that for F ^ imposing the periodic boundary condi- 
tions requires the gauge transformation, which can be 
also seen as interaction representation with respect to 
the Stark term in the Hamiltonian ([2]) . We have checked 
that this procedure does not introduce artifacts and all 
numerical results converge with an increase of the lattice 
size L. 

The solid line in the upper panel in Fig. [T] shows time 
evolution of the mean carrier velocity for the initial condi- 
tions given by the ground state of the system (the Bloch 
wave with zero quasimomentum). It is seen that the 
static force induces Bloch oscillations, which decay to a 
finite V after a transient time ~ I/7. Then, by plotting 
V as the function of F, we reproduce the Esaki-Tsu de- 
pendence (results are not shown). The lower panel in 
Fig. [T] depicts quasimomentum distribution at the end 
of numerical simulations, - a nice correspondence with 
the analytical solution (fTH]) is noticed. We also would 
like to mention the crucial effect of forcing on the system 
coherence, which we characterize by the linear entropy 
c(i) = Tr[/5^]. Without static force the system relax into 
the ground state and c{t) approaches unity. Opposite 
to this, if static field is on, the density matrix tends 
to a diagonal matrix in the momentum representation 



FIG. 2: Decay of coherence (left) and normalized matrix ele- 
ments of the stationary density matrix across the main diag- 
onal (right) for L = 64 (dashed line) and L = 128 (solid line). 
The other parameters are the same as in Fig. 1. 

and almost diagonal matrix in the coordinate represen- 
tation with the diagonal elements pi^i « 1/L. Thus c{t) 
approaches 1/L, as it is illustrated in the left panel in 
Fig.[2l Nevertheless, it would be wrong to say that static 
force makes the system completely incoherent. A weak 
coherence inherited from coherent Bloch oscillations is 
preserved in the form of small but nonzero off-diagonal 
elements pi^m with characteristic profile shown on the 
right panel of Fig. [2] This profile is compared with the 
simplest estimate 

_ 1 1 
^''^ ~ Zl + i(F/7)(Z-m) ' 

which one obtains from ([T]) and ([3]) by setting the hopping 
matrix element in the Hamiltonian ([2|) to zero. 

4- The main drawback of the model considered above 
is that it predicts discontinuous distribution function for 
the carrier velocity (which also holds for a finite temper- 
ature). Next we analyze a more realistic model, which is 
free from this drawback, ft is instructive to begin with 
the case of infinite temperature ^)3\, where the master 
equation can be obtained from the first principles. For 
example, if the bath consists of Bose atoms in an optical 
lattice deep enough to justify the Bose-Hubbard model 
and carriers (for example, Fermi atoms, as in Ref. [1]) in- 
teract with Bose atoms according to Hint = ^X)/ 
the relaxation term in the master equation ([1]) has the 
form j7|], 

i^{p)l,m = - 5Lm)pi,m (14) 

in the Wannier basis, or 

'C(/5)fc,p = -lPk,p + J^}_^ Pk+q,p+q (15) 
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FIG. 3: Numerical analysis of the model No. 2. Stationary 
distributions p{k) for different values of F. Parameters are 
7 = 0.08, J/ksT = 10, and L = 128 (periodic boundary 
conditions) . 
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FIG. 4: Net current for J/ksT = 10, J/fcsT = 1, and 
J/ksT = 0.1 (asterisk), fitted by the Esaki-Tsu equation 
(solid lines). The other parameters are the same as in Fig. 3. 
The inset shows maximally possible current, as compared with 
the temperature dependence f{/3) in Eq. ((Sjl. 



in the Bloch basis, where 7 ~ e^n^ /TlJb with n being the 
density of Bose atoms, Jb hopping matrix elements for 
Bose atoms, and e the interaction constant for collision- 
like interactions of Bose atoms with the carriers [m] . It is 
easy to see that the relaxation operator corresponds 
to the choice W(s,g) = 1/L in the Lindblad form ([6]). 
As a solution of the master equation with relaxation term 
(fT5| one gets the decaying Bloch oscillations and uniform 
stationary distribution for the carrier velocity. 

Since W{s,q) = 1/L implies equal exchange rates 
for any pair of the quasimomentum states, a reasonable 
model for a finite temperature is the set of W^(s, q) where 
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transitions from any lower to any upper energy level are 
suppressed by the Boltzmann factor according to Eq. ([S]). 
Figure [3] shows the stationary distributions p{k) calcu- 
lated on the basis of this model. Unlike in the previously 
considered case W{s,q) — Sg.o, the stationary distribu- 
tions are now smooth functions of k, with the center of 
gravity shifted to positive quasimomenta (F > 0). In ad- 
dition to Fig. [21 Fig. 0] depicts the mean carrier velocity 
as the function of the static force. It is seen that ob- 
tained dependencies are well approximated by the Esaki- 
Tsu formula ([5|), although they do not coincide with ([5|) 
exactly. In particular, with temperature decrease the po- 
sition of the maximum is found to move towards smaller 
F, while in the standard model it is always at F — g. One 
can also see small deviations from the standard model in 
the value of maximal current, shown in inset of Fig. 3) 

5. In conclusion, we have analyzed the Esaki-Tsu equa- 
tion for the drift current with respect to ongoing experi- 
ments on atomic current in the tilted or accelerated op- 
tical lattices. It is shown that different models of carrier 



interactions with the bath result in essentially the same 
(Esaki-Tsu) dependence of the net current on the static 
force magnitude and the bath temperature. This proves 
one more time the universal character of the Esaki-Tsu 
equation. However, different models predict drastically 
different distribution functions for the carriers velocity. 
Since velocity distribution can be directly measured in 
the laboratory experiments with cold atoms, this addi- 
tional information can be used for identifying the details 
of carrier-bath interactions. 

We also mention that the presented analysis questions 
the common believe that the drift current is obligatory 
a non-Markovian process 0, H, UH, where the system 
keep memory about its initial state. Indeed, the above 
discussed carrier dynamics is explicitly Markovian, nev- 
ertheless, it results in the Esaki-Tsu dependence for the 
drift current. 
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